Quantum simulation of the wavefunction to probe frustrated Heisenberg spin systems 
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Quantum simulators are controllable quantum systems that can reproduce the dynamics of the system of in- 
terest, which are unfeasible for classical computers. Recent developments in quantum technology enable the 
precise control of individual quantum particles as required for studying complex quantum systems. Particularly, 
quantum simulators capable of simulating frustrated Heisenberg spin systems provide platforms for understand- 
ing exotic matter such as high- temperature superconductors. Here we report the analog quantum simulation of 
the ground-state wavefunction to probe arbitrary Heisenberg-type interactions among four spin- 1/2 particles . 
Depending on the interaction strength, frustration within the system emerges such that the ground state evolves 
from a localized to a resonating valence-bond state. This spin- 1/2 tetramer is created using the polarization 
states of four photons. The single-particle addressability and tunable measurement-induced interactions provide 
us insights into entanglement dynamics among individual particles. We directly extract ground-state energies 
and pair- wise quantum correlations to observe the monogamy of entanglement. 



During the past years, there has been an explosion of in- 
terest in quantum-enhanced technologies. The applications 
are many-fold and reach from quantum metrology 1 1 ] to quan- 
tum information processing! 2]. In particular quantum com- 
putation has generated a lot of interest due to the discovery 
of quantum algorithms BHH which outperform classical ones. 
The first proposed application for which quantum computa- 
tion can give an exponential enhancement over classical com- 
putation was suggested by Richard Feynman[6, 7]. He con- 
sidered a universal quantum mechanical simulator, which is a 
controllable quantum system that can be used to imitate other 
quantum systems, therefore being able to tackle problems that 
are intractable on classical computers. Since then the moti- 
vation to use a quantum simulator as a powerful tool to ad- 
dress the most important and difficult problems in multidisci- 
plinary science has led to many theoretical proposals 18HT31. 
Vast technological developments allowed for recent realiza- 
tions of such devices in atoms fT4lfT6lL trapped ions fT7lj20lL 
single photons EHS) and NMR ll25ll26l . The quantum simu- 
lation of strongly correlated quantum systems (e.g. frustrated 
spin systems) is of special interest and would provide new re- 
sults that cannot be otherwise classically simulated! 27 ]. 

In order to manipulate and measure individual properties of 
microscopic quantum systems the complete control over all 
degrees of freedom for each particle is required. Typically, 
atoms in optical lattices! 14] are used for realizing physical 
systems that can simulate various models in condensed-matter 
physics. The fact that the experimental addressability of sin- 
gle atoms in optical lattices remains very challenging l28ti30l 
leads to the studies of bulk properties of the atomic ensemble 
10 5 atoms) instead of single particles. Therefore we utilize 
single photons in separate spatial modes and measurement- 
induced interactions as a quantum simulator, thus the parti- 
cles are individually accessible. The tunable interaction be- 



tween two entangled photon-pairs allows for the precise sim- 
ulation of the ground state of a spin- 1/2 tetramer. We obtain 
the ground- state energy and have direct access to the distri- 
bution of pair- wise quantum correlations as a function of the 
competing spin- spin interactions. We also observe the influ- 
ence of mono gamy OTl l32l in this strongly correlated quantum 
system. 

Analog quantum simulator 

The main challenge in the understanding of strongly corre- 
lated quantum systems is to calculate the energies and ground 
state properties of many-body systems as this becomes expo- 
nentially difficult with increasing number of particles when 
using a classical computer. In contrast, quantum simulators 
use quantum systems to store and process data which allows 
them to polynomially mimic the evolution of the quantum 
system of interest. Typically, quantum simulations require 
methods to implement the Hamiltonian of the simulated sys- 
tem, and probe its ground state properties. In some cases, 
even if the ground state wave-function is known, obtaining all 
properties by using classical computers remains a challeng- 
ing task. For example, for the one-dimensional antiferromag- 
netic Heisenberg model one can compute the spectrum and 
the wave-function by using the Bethe ansatz l33lL but it is not 
as easy to extract the correlations functions. This can be over- 
come by quantum simulators that are capable of generating 
and directly probing the quantum state of interest via control- 
ling each quantum particle individually. 

Usually, the system being simulated is defined by its Hamil- 
tonian H(t,J,B,...) that is dependent on parameters such as 
time, t, interaction strength, J, external field, B, etc. In prac- 
tise, one way of realizing a quantum simulator is based on dis- 
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FIG. 1 : Mimicking an adiabatic quantum evolution with an analog 
quantum simulator, (a) Adiabatic quantum evolution. The system 
is prepared in an initial ground state \ip(to, Bo, Jo, . . . )). Then the 
gradual change of the system parameters (£, B, J, etc.) causes an 
adiabatic evolution of the system to the final ground state of inter- 
est \ip(t, B, J, ...))• ( D ) Analog quantum simulation. The adiabatic 
evolution of the system to be simulated is mapped onto a controllable 
evolution of a quantum system. A set of tunable gates give access to 
the change of parameters, (c) Model used to study the valence-bond 
states. The nearest-neighbor Heisenberg-type interactions of strength 
Ji and J 2 among four spin- 1/2 particles are drawn as connecting 
bonds and form a spin- 1/2 tetramer. All the properties of the tetramer 
depend only upon the coupling ratio k = J2/J1. (d) Quantum simu- 
lation of a spin- 1/2 tetramer using a photonic analog quantum simu- 
lator. The initial ground state, |\P(0o)), is prepared by generating the 
photon-pairs 1 & 2 and 3 & 4 in two singlet states. Then the analog 
quantum simulation is performed utilizing the measurement-induced 
interaction, consisting of quantum interference and the detection of 
a four-photon coincidence after superimposing photons 1 & 3 on a 
tunable beam splitter. Mapping the coupling ratio k on the beam 
splitter's splitting ratio tan 2 leads to the ground state of interest, 
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crete gate operations and the phase estimation algorithmfTOl 
[231 , referred to as a digital quantum simulator! 12]. An alter- 
native approach utilizes the adiabatic theorem! 34], where an 
initial Hamiltonian, whose ground state is easy to prepare, can 
be adiabatically evolved to a final Hamiltonian with a nontriv- 
ial ground state of interestlHl l9l 1131 126) . An adiabatic quan- 
tum simulator can be built by engineering interactions among 
particles using tunable external parameters (e.g. an external 
magnetic field). The system will remain in its ground state if 
the system parameters change gradually enough. 

Our experimental technique combines the advantages of 
both approaches by utilizing a tunable quantum gate without 
the necessity of either discretizing the quantum evolution or 
engineering the physical interactions for an adiabatic quan- 
tum simulation. Thus, we consider our simulator as an analog 
quantum computer |[T2l [T4l [T8l l20l [28ti30l , where the change 
of the quantum evolution can be obtained by a tunable quan- 
tum gate. The advantages of using analog quantum gate op- 
erations seem to be far-reaching as the number of gate oper- 
ations can be less-demanding than for digital quantum simu- 
lators. Fig. [T] shows the concept of this analog simulator. In 
our experiment, we use the analog quantum simulator to di- 



rect probe the wave-function instead of the Hamiltonian. We 
take the advantage of single-particle addressability and the 
tunable interaction of our quantum simulator to prepare vari- 
ous ground states experimentally. Once the wave-function is 
prepared, it is feasible to measure all the correlation functions 
and study the dynamics of entanglement in this strongly cor- 
related systems. 



Simulation of a spin- 1/2 tetramer 

Over the last 60 years the fundamental interest in stud- 
ies of ground states of Heisenberg-type Hamiltonians has 
led to a few exact theorems, which may serve as guidelines 
for quantum simulators. Based on Marshall's theorem l35l 
and its extension l36l the absolute ground state has total 
spin zero (S 2 = 0) for N spins on a bipartite lattice with 
nearest-neighbor antiferromagnetic Heisenberg-type interac- 
tions. This constraint leads to the fact that the ground state can 
be built as a superposition of pairs of spins forming singlets or 
equally so-called valence bonds. If all the spins are covered 
by valence bonds, which are maximally entangled states, then 
the ground state's total spin is zero and non-magnetic. This 
is established by valence bonds that are either static and lo- 
calized or fluctuating as a superposition of different partition- 
ings of spins. In general, the equally weighted superposition 
of all different localized valence-bond states corresponds to 
a quantum spin liquid, the so-called resonating valence-bond 
state|f37][38l. 

The smallest configuration for studying and simulating 
these phenomena on a two-dimensional square lattice is four 
spin- 1/2 particles forming a tetramer. In the case of such a 
spin- 1/2 tetramer the Heisenberg-type interactions lead to the 
creation of three possible dimer-covering configurations for 
the localized valence-bond states, |<£=) = |^ _ ) 12 |^~) 34 , 

\®\\) = l^~>13^~>24 and l$x> = I^~ii4 l^~>23» Where 

is me singlet of particle i and j (Fig. fit). Since the total 
spin-zero subspace for this system is two-aimensional, these 
three dimer-covering states are not independent and |3> x ) 
can be written as |$ x ) = |$y) - |<£=) in the |4>=)/|$||) 
basis. This state, |$ x ), like any other equal superposition 
of these two dimer-covering states represents a resonating 
valence-bond state. Particularly interesting states are resonat- 
ing valence-bond states! 1 1 . 40] with s-wave pairing symmetry, 
I $ || ) + 1 3>= ) (up to normalization), and with the exotic d-wave 
pairing symmetry, |$ x ) = |$||) — |$ = ). The studies of these 
states are of high interest, because it was conjectured that a 
transition from an localized valence-bond configuration to the 
superposition of different valence-bond states might explain 
high-temperature superconductivity in cuprates|39|. A quan- 
tum simulator capable of preparing such arbitrary superposi- 
tions of dimer-covering states is thus sufficient for simulating 
any Heisenberg-type interactions of four spin- 1/2 particles on 
a two-dimensional lattice. It is the particular strength of our 
optical quantum simulator that the simulated ground states can 
be restricted to the spin-zero singlet subspace by utilizing the 



quantum interference of photons at a tunable beam splitter. 

For our quantum simulation we consider four spin- 1/2 
particles on a square lattice (tetramer) that interact via 
Heisenberg-type interactions. According to Marshall's the- 
orem, the ground state is restricted to total spin-zero subspace 
for arbitrary antiferromagnetic interactions, independent of 
their type and strength. Here we model our spin tetramer with 
nearest-neighbor interactions of the strength J\ and J 2 (Fig. 
[ij) by the Hamiltonian 

H = J1S1S3 + J1S2S4 + J2S1S2 + J2S3S4, 

where Si is the Pauli spin operator for spin i. All the prop- 
erties of the system depend only on the coupling ratio k = 
J 2 1 J\ therefore we re-normalize the Hamiltonian to 

H{k) = H + kH 1 , (1) 

where H(k) = Hj J\ is the final Hamiltonian, Hq = S1S3 + 
S2S4 the initial Hamiltonian and Hi = S1S2 + S3S4 the 
competing Hamiltonian of Ho. 

Due to the simplicity of the quantum system we study here, 
we can calculate the expected ground state analytically, which 
allows us to verify the experimental results. The ground state 
of the Hamiltonian given in Eq. ([T]) is 

*<°>(0)) = _l-(cos2fl|$ = > -cos 2 ^,)) 

1 V n (°) 

where n{0) = ^(cos 4 + cos 2 20 + sin 4 0) is the normaliza- 
tion constantThe ground state energy of \^^\0)) is = 
-2(1 + k) - 4Vl-^ + ^ 2 . The TDC's angle, 0, takes the 
values from the interval of < < | . Using our photonic 
quantum simulator, we can mimic the adiabatic change of the 
Hamiltonian shown in Eq. ([T]), where the full range of the cou- 
pling ratio —00 < k < +00 is experimentally covered by 
tuning the angle of the TDC between arctan < < | . 

For k = 0(0 = f ), the ground state is \®\\), while for 
k = +00 (6 = |), the ground state changes to |3>=). These 
two cases are dimer-covering states. Tuning the coupling ra- 
tio to k = —00 results in the equally weighted superposition 
|3>x) + I $|| )» whereas k = 1 leads to the interesting resonat- 
ing valence-bond state t37ti39l |$ = ) + |$y). 

Here, we experimentally demonstrate an optical analog 
quantum simulator by producing two polarization entangled 
photon pairs (see Fig. go, \^~) 12 = ^(|HV) 12 - |VH) 12 ) 

and = ^(|HV) 34 — |VH) 34 ), in the spatial modes 1 

& 2 and 3 & 4. \H) and | V) denote horizontal and vertical po- 
larization states, respectively. The tunable interaction among 
these singlet states is achieved by the TDC, followed by a pro- 
jective measurement of one photon in each of the four output 
modes. This tunability allows us to continuously change the 
measurement-induced interaction between photons 1 and 3. 
The TDC is an optical fiber device that transfers optical sig- 
nals between fibers acting as a beam splitter with controllable 
splitting ratio. The control of the splitting ratio is achieved 
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FIG. 2: Experimental setup, (a) Femtosecond laser pulses 140 fs, 
76 MHz, 404 nm) penetrate two /3-barium borate (BBO) crystals 
generating two pairs of photons in the spatial modes 1 & 2 and 3 & 
4 (two-fold coincident count rate per pair w 20 kHz). The walk-off 
effects are compensated with a half- wave plate (HWP) followed by a 
BBO crystal in each mode. The photon's spectral and spatial distin- 
guishability is erased with interference filters (IF, FWHM = 3 nm) 
and single-mode fibers. The polarization of each photon is analyzed 
by a combination of a quarter- wave plate (QWP), a HWP and a po- 
larizing beam splitter (PBS). Single photons are detected by single- 
photon counting modules (SPCM). (b) Schematic diagram of the 
fiber-based tunable directional coupler (TDC). The view from the top 
of the TDC illustrates the dependence of the coupling of the evanes- 
cent light upon the separation of the fibers. The coupling between 
these two fibers is controlled by adjusting the horizontal position of 
the D fiber, (c) Experimental calibration of TDC's transmissivity (red 
circles) and reflectivity (black circles) with respect to the position of 
the D fiber (s) is performed by using weak laser beams and SPCM. 
The separations of the fibers for 0%, 50% and 100% transmissivity 
are shown in the insets. The experimental imperfections originate 
mostly from detector dark counts and the error bars, smaller than 
0.5% of the mean values, are based on a Poissonian distribution. 



by adjusting the relative positions of the fibers (Fig. [2J3). The 
transmissivity and reflectivity of this TDC go as cos 2 and 
sin 2 0, respectively, where parameterizes the separation of 
the fibers. We calibrate the TDC's transmissivity and reflectiv- 
ity such that the modulating visibilities (Michaelson visibility) 
are above 95% for both inputs, as required for high-precision 
quantum control (Fig. [2]:). 

A successful detection of a four-fold coincidence event 
from each spatial mode gives the four-photon state, 

\m)i2M = 7fe[- cos2 ^i HHVV > + |WHH» + 

sin 2 0(|HWH) + |VHHV)) +cos 2<9(|HVHV) + |VHVH))]. 
The experimentally obtained density matrix, p exp , is recon- 
structed from a set of 1,296 local measurements using the 
maximum-likelihood technique BTl l42l . For this, all combi- 
nations of mutually unbiased basis sets for individual qubits, 
that is |H/V>, |+/-) = ^(|H> ± |V» and |R/L> = 

^(|H) ± i |V)), are measured. The duration of each mea- 
surement for a given setting of the polarization analyzers and 
the TDC was 200 s and the average detected four-fold coin- 
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FIG. 3: Ground state energy of the spin- 1/2 tetramer. By tuning 6, 
where tan 2 — k + V k 2 — k + 1 represents the splitting ratio of 
the tunable directional coupler, we gradually change the ground state 
of the spin- 1/2 tetramer. The full range of the coupling ratio — oo < 
k — < oo is covered by tuning from arctan 4^ to \ . We 
measure the ground state energy for seven different configurations. 
Of particular interest are the quantum states |$y) + |$x)» |^||)» 
| $H ) + |3>=) and |$=), shown explicitly. The black circles represent 
the experimental data and the solid line is parameter-free theoretical 
prediction. The error bars follow Poissonian statistics and are smaller 
than the data points. 



FIG. 4: Density matrices of various spin- 1/2 tetramer configura- 
tions in the computational basis (\H)/\V)). Shown are the real parts 
of the density matrices for the cases of (a) equal superposition of 
dimer-covering states, (b,d) dimer-covering states and (c) resonat- 
ing valence-bond state. The imaginary parts are small and shown in 
the Supplementary Information. The wire grids indicate the expected 
values for the ideal case. The density matrices are reconstructed from 
the experimental four-photon tomography data for the settings of (a) 
= 0.197tt, (b) = 0.25tt, (c) = 0.304tt and (d) = 0.468tt. 
The fidelities, F, of the measured density matrix with the ideal state 
are (a) F = 0.745(4), (b) F = 0.712(4), (c) F = 0.746(6) and (d) 
F — 0.888(2). The uncertainties in fidelities extracted from these 
density matrices are calculated using a Monte Carlo routine and as- 
sumed Poissonian errors. 



cidence rate was 3 Hz. In total eight density matrices for dif- 
ferent settings of 6 are reconstructed and are summarized in 
the Supplementary Information. Uncertainties in quantities 
extracted from these density matrices are calculated using a 
10 run Monte Carlo simulation of the whole state tomography 
analysis, with Poissonian noise added to each experimental 
data point in each run. 

The ground-state energy, E^°\ is defined as the mean value 
of the Hamiltonian given in Eq. ([I]). Since this Hamiltonian 
contains only pair- wise interactions, we measured the expec- 
tation value of the corresponding pair-wise correlations and 
obtained the data shown in Fig. |3j which show good agree- 
ment with theoretical prediction. 

In Fig|4j we show the experimentally obtained density ma- 
trices of the four valence-bond states, |3>=) + |3> x ), 3>||), 
|3>=) + \<&\\) and |3>=), which correspond to the setting of 
= 0.197tt, = 0.25tt, = 0.304tt and = 0.468tt. The 
state fidelity is defined as F(*,p) = (*|p|*), where \Sfr) 
is the target and p is experimentally obtained quantum state. 
Due to the high quality of our quantum simulator, we obtain 
four-photon state fidelities that range from F = 0.712(4) to 
F = 0.888(2) (see Supplementary Information). 



Quantum monogamy and complementarity 

Monogamy is one of the most fundamental properties of 
quantum entanglement | 3TJ [32). It restricts the shareability 
of quantum correlations among parties and is of essential im- 



portance in many quantum information processing protocols, 
including quantum cryptography and entanglement distilla- 
tion. Recent work showed that in the context of condensed- 
matter physics, monogamy gives rise to frustration effects in 
quantum system e.g. Heisenberg antiferromagnets. The ideal 
ground state for an antiferromagnet would consist of singlets 
between all interacting spins. But, due to the monogamy re- 
lation a particle can only share one unit of entanglement (sin- 
glet) with its neighbors. Therefore, it will spread entangle- 
ment in an optimal way with all its neighbors leading to a 
strongly correlated ground state (32). 

To study the dynamics of pair- wise interactions, in which 
the monogamy of bipartite quantum entanglement distribu- 
tion plays a crucial role, we characterize the distribution of the 
two-body energies and correlations between one spin with re- 
spect to the others with the normalized Heisenberg energy per 
unit of interaction, e^ . It is defined as = — ^Tr(pijSiSj). 
Note that pij is the density matrix of spins i and j. The 
normalized Heisenberg energy per unit of interaction is also 
an entanglement witness 1 43, 44] and reaches its maximum 
value of eij = 1 for the singlet state. The amount of en- 
tanglement can also be quantified by concurrence (45), which 
is directly related to with C(e^) = max{0, — \ + |e^}. 
For our four- spin system the dependencies of the pair- wise 
energies with respect to the TDC's angle are given by 
ei2 = — ^(sin 2 6 ) cos26 ) ), ei3 = ^(sin 2 cos 2 0), and eu = 
^(cos 2 <9cos2<9). 

Remarkably, monogamy is manifested in the constraint of 
the energy distribution for the considered spin pair through a 
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FIG. 5: Experimentally extracted pair- wise Heisenberg energies, (a) 
Experimental observation of quantum monogamy when comparing 
the pair- wise normalized Heisenberg energy, e^ , it acts as a two- 
particle entanglement witness, for the spin pairs 1 & 2 (black square), 

1 & 3 (red circle) and 1 & 4 (blue triangle). The highlighted area cor- 
responds to the full range of the coupling ratio oo. 
For the case of k = (0 — J), the ground state of this spin- 1/2 
tetrameris \&\\) — |^ _ ) 13 \^~) 2 4 an d the amount of entanglement 
of the pair 1 & 3 reaches its maximum while the pairs 1 & 2 and 1 & 
4 are not entangled. Similarly, for the case of k — +oc (0 = f), the 

ground state is reduced to |<3>=) = \^~) 12 \^~) ^ wnere P a i r 1 & 

2 is now maximally entangled, and pairs 1 & 3 and 1 & 4 are disen- 
tangled. In the case of the resonating valence-bond state, entangle- 
ment distributions are equal between the pairs 1 & 2 and 1 & 3 (i.e. 
ei2 = eis). In other cases, entanglement is distributed according 
to the monogamy relation, (b) Experimental demonstration of the 
complementarity relation in a spin- 1/2 tetramer. For each valence- 
bond configuration we measured pair- wise Heisenberg energies, e^ , 
which are normalized by its maximal value, e™ ax . The sum of these 
renormalized energy values are in good agreement with the theoreti- 
cal prediction (shown as line in the plot). The uncertainties represent 
standard deviations deduced from propagated Poissonian statistics. 
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This restricts the maximal amount of energy or entanglement 
associated with correlated spin systems (see Fig. [5^). For in- 



stance, in the experiment we obtain the normalized Heisen- 
berg energy per unit interaction between photons 1 and 2(eyi) 
with the correlation measurements in three mutually unbiased 
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FIG. 6: Directly observed pair- wise correlation functions of various 
valence-bond states. The correlation tensors T\ 2 (photons 1 & 2), 
T13 (photons 1 & 3) and T3.4 (photons 1 & 4) are obtained from cor- 
relation measurements directly in the bases X = a x ,Y = a y and 
Z — a z . For a convenient graphical representation, the negative 
values of the correlation tensors are shown. The structure of (a) the 
superposition state and (c) resonating valence-bond state show that 
the quantum correlations are equally distributed among two compet- 
ing pairs, (b) and (d) belong to dimer-covering states, in which only 
one pair is maximally correlated in a singlet state. 
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As shown in Fig. [5^, the adiabatic change of the coupling 
between the four spins is simulated by tuning the angle of the 
TDC, 0, from arctan ^= to | . This corresponds to the full 

range of the coupling ratio —00 < k = ^ < 00. In the 
ideal case, the maximum of is unity which corresponds to 
a singlet state shared by spins i and j. However, imperfec- 
tions in the generation of entangled photon pairs and the two- 
photon interference on the TDC reduce the measured value of 
eij by a constant factor, independent of 0. For the individ- 
ual photon pairs we obtain the maximal Heisenberg energy of 
e£f x = 0.920(7), e™ 3 ax = 0.727(9), and ef£* = 0.926(5). 
In order to demonstrate the complementarity relation (48) we 
re-normalized each energy by its maximal value e™ ax and 
obtain a good agreement with the theoretical prediction shown 
in Fig(5j3. 

The advantage of the individual addressability for our par- 
ticles in the ground state allows for the direct extraction of 
the pair- wise quantum correlations. The pair- wise quantum 
correlation is defined as: 
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where U ) are the coincidence counts between 

pair i and j along the direction of and that perpendicular 
to S( v \ respectively. In Fig.[6j the pair- wise correlation func- 
tions for the ground states, |$ = ) + |$x), |$||), |$=) + |^||) 

and |3>=), are shown. As expected from the monogamy rela- 
tion, in the cases of dimer-covering states, one pair of the pho- 
tons is maximally correlated, e.g. photons 1 and 3 in Fig. [6]), 
and photons 1 and 2 in Fig. |6j:. In the cases of the equal super- 
position of two dimers (Fig. [6^) and the resonating valence- 
bond state (Fig.[6]l), correlations are distributed among differ- 
ent pairs. 

We demonstrate the feasibility of an all-optical analog 
quantum simulator by enabling quantum control of the 
measurement-induced interaction among photonic quantum 
states. Various ground states, including the resonating 
valence-bond states for four interacting spin- 1/2 particles are 
generated and characterized by extracting the total energy 
and the pair- wise quantum correlations. The simulation of a 
spin- 1/2 tetramer also proves that the pair- wise entanglement 
and energy distribution are restricted by the role of quantum 
monogamy. Our results provide promising insights for quan- 
tum simulations of small quantum systems, where individual 
addressability and control over all degrees of freedom on the 
single-particle level is required. This is of particular inter- 
est for quantum chemistry with small numbers of particles 
and might allow in the near future the simulation of aromatic 
systems and chemical reactions (49). Although it was shown 
that efficient scalable quantum computing with single pho- 
tons, linear-optical elements, and projective measurements is 
possible l50lL the most important challenges will be (a) the 
realization of two- and more-qubit interactions with high fi- 
delity, (b) generating systems with more qubits and (c) devel- 
oping efficient methods of simulating other classes of com- 
plex Hamiltonians by using optical elements. Ideally, this and 
related work will open a new and promising avenue for the 
experimental simulation of various quantum systems. 
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248095), an ERC senior grant (QIT4QAD), the Marie-Curie 
research training network EMALI, JTF, SFB-FOQUS and the 
doctoral programme CoQuS of the Austrian Science Founda- 
tion (FWF). 
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Supplementary Information 
Quantum state tomography 

We perform quantum state tomography for the simulated 
ground states of a spin- 1/2 tetramer. The studies of the eight 
different ground states lead to 10368 coincidence count mea- 
surements with a total 414724 four-fold coincidence counts. 
The reconstructed density matrices are plotted in Fig.[7]and[8] 





FIG. 8: Density matrices of various spin- 1/2 tetramer configurations 
in the computational basis (\H)/\V)). Shown are the real parts (p re ), 
imaginary part (p l7n ), and absolute values (p abs ) of the density ma- 
trices for the different settings of the splitting ratio of of the tunable 
direction coupler: (a), 6 = 0.25tt, (b), 6 = 0.222tt, (c), 6 = 0.197tt 
and (d), = 0.0457T. The wire grids indicate the expected values 
for the ideal case. The fidelity, F, of the measured density matrix 
with the ideal state are (a), F = 0.712(4), (b), F = 0.755(5), (c), 
F = 0.745(4) and (d), F = 0.897(1). 



FIG. 7: Density matrices of various spin- 1/2 tetramer configurations 
in the computational basis (\H)/\V)). Shown are the real parts (p re ), 
imaginary part (p l7n ), and absolute values (p abs ) of the density ma- 
trices for the different settings of the splitting ratio of of the tun- 
able direction coupler: (a), 9 = 0.4687T, (b), = 0.4557T, (c), 
= 0.3667T and (d), = 0.3047T. The wire grids indicate the 
expected values for the ideal case. The fidelities, F, of the mea- 
sured density matrix with the ideal state are (a), F = 0.888(2), (b), 
F = 0.888(2), (c), F = 0.840(4) and (d), F = 0.746(6). 
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